Failure of Mean Field Theory at Large N 
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We study strongly coupled lattice QCD with TV colors of staggered fermions in 3 + 1 dimensions. 
While mean field theory describes the low temperature behavior of this theory at large N, it fails in 
the scaling region close to the finite temperature second order chiral phase transition. The universal 
critical region close to the phase transition belongs to the 3d XY universality class even when TV 
becomes large. This is in contrast to Gross-Neveu models where the critical region shrinks as TV 
(the number of flavors) increases and mean field theory is expected to describe the phase transition 
exactly in the limit of infinite TV. Our work demonstrates that close to second order phase transitions 
infrared fluctuations can sometimes be important even when TV is strictly infinite. 



Mean field techniques provide a simple but powerful 
approach to gain qualitative insight of the underlying 
physics in a variety of field theories The Bardeen- 
Cooper-Schrieffer (BCS) solution to superconductivity 
is a well known application of such a technique. Wil- 
son's renormalization group shows that when correlation 
lengths £, associated with the fluctuations of the field, 
become large compared to the microscopic length scale 
a of the problem, mean field techniques become exact 
in dimensions greater than four. However, in lower di- 
mensional systems infrared fluctuations can become im- 
portant when £ 3> a and invalidate the mean field ar- 
guments. For this reason the mean field approach must 
be used with care close to second order phase transi- 
tions. The region close to the critical point where mean 
field theory fails is usually referred to as the Ginsburg 
region Q. Using field theoretic techniques sometimes it 
is possible to estimate the Ginsburg region. In conven- 
tional superconductors the Ginsburg region is known to 
be suppressed by some power of T c /Ef where T c is the 
critical temperature and Ef is the Fermi energy. Mean 
field theory is believed to be reliable outside the Ginsburg 
region. This is the reason why BCS is a good approach to 
understand the physics of superconductivity for all tem- 
peratures except very close to T c . 

There are certain limits in which the infrared fluctua- 
tions can be naturally suppressed even in low dimensional 
systems. For example consider a theory containing a field 
with TV components. When TV is large a saddle point 
approximation can be used to solve the theory and the 
leading term is nothing but the mean field solution. The- 
oretical physicists often use this feature to solve a phys- 
ical theory by increasing the number of components of 
the field artificially. By solving the theory in the limit of 
large number of components and computing the leading 
corrections they can sometimes estimate realistic answers 
in the physical theory. A variety of field theories can be 
studied in the large TV limit [3j- 

Although large TV approach to field theories has been 
very successful, in this article we show that not all large 



TV limits lead to the suppression of infrared fluctuations. 
We contrast two models, the Gross-Neveu (GN) model 
involving TV species of fermions studied recently 0, Q 
and strongly coupled U(N) lattice QCD with staggered 
fermions (SCLQCD) studied here. Both these models 
contain a global symmetry which is spontaneously bro- 
ken at low temperatures. Further, the symmetry group 
and the breaking pattern is not affected by TV. These 
models undergo a second order phase transition to a sym- 
metric phase at a finite critical temperature T c . In the 
large TV limit they can be solved exactly using mean field 
techniques at low temperatures. An interesting question 
then is whether the mean field description is valid even 
close to T c when TV becomes infinite. 

It was discovered in 0] that indeed in the GN model 
the critical behavior near T c belongs to the Landau- 
Ginsburg mean field universality class at large TV. This 
implies that the Ginsburg region, where the critical be- 
havior is governed by a non-trivial universality class that 
depends on the symmetry group and the breaking pat- 
tern in a dimensionally reduced theory, has a zero width 
at large TV. It was later shown in [j| that indeed the Gins- 
burg region is suppressed by a factor Although 
a Z2 symmetric GN model was analyzed in 0, [E[ , there 
are reasons to believe that the arguments would hold for 
continuous symmetries as well jfj, |j| • 

As we will demonstrate here, in contrast to the GN 
model, the Ginsburg region does not shrink in SCLQCD 
in the large TV limit. While mean field theory is indeed a 
good approximation at low temperatures, the finite tem- 
perature phase transition is not described by mean field 
theory even at infinite TV. We believe our results should 
be of interest to a wide range of physicists since our model 
can be mapped into a theory of classical dimers. Dimcr 
models have a long history In the 1960s these models 
attracted a lot of attention when it was shown that the 
Ising model can be rewritten as a dimer model 
the late 1980s they gained popularity again in their quan- 
tum version \\\\ as a promising approach to the famous 
resonating- valence-bond (RVB) liquid phase More 
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recently, this approach has gained momentum again since 
it was shown that the RVB phase was indeed realized on 
a triangular lattice but not on a cubic lattice 0, Hif . 
Thus physicists attempting to use large N techniques in 
dimer models can benefit from our results. 

The partition function of SCLQCD is given by 

Z(T,m} = J \dU][dipdip] exp (S[U, ip, ip]) , (1) 

where [dU] is the Haar measure over U(N) matrices and 
[dip dip] specify Grassmann integration. The Euclidean 
space action S\U,ip,ip] in the strong (gauge) coupling 
limit with staggered fermions is given by 
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where x refers to the lattice site on a periodic four- 
dimensional hyper-cubic lattice of size L along the three 
spatial directions and size Lt along the euclidean time di- 
rection. The index /i = 1, 2, 3, 4 refers to the four space- 
time directions, U Xifi G U(N) is the usual links matrix 
representing the gauge fields, and ip x ,ipx are the three- 
component staggered quark fields. The gauge fields sat- 
isfy periodic boundary conditions while the quark fields 
satisfy either periodic or anti-periodic boundary condi- 
tions. The factors rj x ^ are the well-known staggered 
fermion phase factors. Using an asymmetry factor be- 
tween space and time we introduce a temperature in the 
theory. We choose r] x „ = 1, fi — 1,2,3 (spatial direc- 
tions) and r] 2 4 = T (temporal direction), where the real 
parameter T is a coupling that controls the tempera- 
ture. By working on anisotropic lattices with L t << L at 
fixed L t and varying T continuously one can study finite 
temperature phase transitions jl. r | . In this article we fix 
L t = 4 for convenience. 

The partition function given in eq.0 can be rewrit- 
ten as a partition function for a monomer-dimer system, 
which is given by 
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and is discussed in detail in 0, [l^. Here n x = 
0, 1, 2, .., N refers to the number of monomers on the site 
x, b x _ fj, — Q, 1,2, ...,N represents the number of dimers 
on the bond connecting x and x + fi, to is the monomer 
weight, z x ,fi = ?7x,u/4 are the dimer weights. Note that 
while spatial dimers carry a weight 1/4, temporal dimers 
carry a weight T/4. The sum is over all monomer-dimer 
configurations [n, b] which are constrained such that at 
each site, n x + J^^fix^ + b x -p,,n] = N. 

When to = 0, the action of SCLQCD, cq. QJ, 
is invariant under 0(2) chiral transformations: ip x — > 
exp (ia x 0)ip x and ip x — > ip x eicp(ia x 9) where a x = 1 for all 
even sites and —1 for all odd sites. In the large N limit 



mean field techniques can be used to show that this chiral 
symmetry is spontaneously broken at low temperatures 
flfij |. In a Monte Carlo method was developed to 
solve the problem from first principles and it was shown 
that mean field theory is indeed reliable at small tem- 
peratures 01- Unfortunately, since the algorithm was 
inefficient at small quark masses, the finite temperature 
chiral phase transition was never studied in the large N 
limit. Recently a very efficient cluster algorithm was dis- 
covered to solve the model at any value of N . Using 
this algorithm it was shown with great precision that for 
TV = 3 the finite temperature phase transition belonged 
to the 3d XY universality class • Here we extend that 
calculation to higher values of TV. 

The order parameter that signals chiral symmetry 
breaking is the chiral condensate, defined by 



lim lim 

m— »0 L—>cg 



1 



1 d 



NL t L 3 Z dm 



Z(T,m) 



(4) 



For a fixed T the large N result for S can easily be ob- 
tained by extending the calculation of [r^. One gets 
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which shows that the critical temperature, as we have 
defined it, is infinite. A calculation which includes the 
1/N correction shows that T c ~ N. In the large d (spatial 
dimensions) limit one obtains T c = d(N + 2)L t /6 |2(l |. 

In order to determine S using Monte Carlo calculations 
we measure the chiral susceptibility in the chiral limit, 
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The finite size scaling of this quantity is known from chi- 
ral perturbation theory |2l| and one expects 
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The constant F 2 is equal to 
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at to = 0. The current J x „ is the conserved current asso- 
ciated with the 0(2) chiral symmetry ^^|. By fitting the 
data for x to the form given in eq. Q we can determine 
S accurately. 

We have done extensive simulations for various values 
of L and N in order to extract E as discussed above. In 
fig. © we plot S as a function of T at N = 6, 12, 24, 48. 
For comparison we also plot the mean field result (eq.jsj)). 
As the graph shows, at a fixed value of T, our data ap- 
proaches the mean field prediction quite nicely as A^ be- 
comes large. However, for every value oi N, as T in- 
creases the order parameter approaches zero at some crit- 
ical temperature T c . Close to T c , the mean field theory 
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where one holds T/N fixed as N becomes large. Let us 
refer to this as the finite T mean field theory (FTMFT). 
Although we have not yet developed this mean field the- 
ory, we will argue that it is bound to fail close to T c since 
the Ginsburg region where the 3d XY universality class 
is observed does not shrink with N. Indeed we find that 
E ~ (1 - T/Trf dose to T c where (3 = 0.3485(2) inde- 
pendent of N [22| . A fit of the data to this form was used 
to determine T c plotted inQ In fig- 0we plot E as a func- 
tion of T/T c for 0.95 < T/T c < 1 for N = 6, 12,24,48. 
The solid lines are fits to the form A(l - T/T c ) - 3485 . 
We find A = 0.3668(3), 0.2630(2), 0.1862(2), 0.1315(2) for 
N = 6,12,24,48 respectively all with a \ 2 /DOF less 
than 1. 



FIG. 1: Plot of E vs. T for various values of N. The errors 
in the Monte Carlo data are less than the size of the symbols. 
The solid line is the mean field result given in eq. (Jjj|- 



is definitely not a good approximation. Using our algo- 
rithm we can determine T c accurately for every value of 
N (see below). In the inset of fig. we plot T c as a 
function of N. A fit to the form T c = aN + b + c/N 
yields a = 1.5525(3), b = 3.126(9) and c = -0.88(4) with 
a x 2 /DOF = 0.73 (solid line in the inset). 




FIG. 2: Plot of E vs. T/T c for various values of N. The solid 
lines represent the fit to A(l - T/T c ) 3485 which is expected 
from 3d XY model. The values of A for different N are given 
in the text. 

One might think that it is quite easy to explain why 
the mean field theory breaks down close to T c . Since T c 
grows as N, the large N theory does not know about 
the existence of a finite T c unless the 1/7V corrections are 
included. We have computed these corrections and have 
found that they do not improve the situation much. Per- 
haps one needs to develop a new mean field expansion 




FIG. 3: Plot of F 2 vs. T/T c for N = 48, 24, 12. The inset 
shows the same graph for a larger range of T/T c for N — 48. 
The text describes the physics of this plot. 

Since A decreases with increasing N one might argue 
that A will vanish in the limit N — ► oo. In that case some 
higher order term will become dominant in the large N 
limit and it is not possible to rule out mean field be- 
havior. However, we do not think this is the case and 
attribute the change in A to a renormalization effect as 
a function of N. In order to justify this, we next focus 
on a correlation length scale which does not need renor- 
malization. If correlation lengths in the theory indeed 
scale with N, we would expect the correlation lengths to 
be the same at any fixed T/T c . In our model F 2 can be 
defined as one such inverse correlation length scale. This 
implies F 2 - B(l - T/T c ) 1 ', where v = 0.6175(3) in the 
3d XY model |22j. In fig. [3] we plot F 2 as a function 
of T/T c for N = 12, 24, 48. When T/T c > 0.8 all the 
points at a fixed T/T c but at different N fall on top of 
each other. Further, when 0.95 < T/T c < 1, all the data 
points fit extremely well to the form B{\ - T/T c ) - 6715 , 
with B = 0.939 with a x 2 /DOF = 0.65. We believe that 
this is strong evidence that even at infinite N the phase 
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transition belongs to the 3d XY universality class. Out- 
side the Ginsburg region one would expect FTMFT to be 
valid 0. Interestingly, when 0.85 < T/T c < 0.93 we find 
that the data fits to the form B(l-T/T c ) with B = 1.507 
with a x 2 /DOF — 0.8, suggesting v = 1. Hence, we sus- 
pect that the FTMFT in our model is similar to the mean 
field theory in 3d O(N) models which yields v = 1. The 
cross over to mean field theory occurs when T/T c < 0.95 
so that the correlation length is £ ~ 1/F 2 < 8 in lattice 
units. 

One often hears the lore that in the large N limit 
SCLQCD is solvable using mean field theory. While this 
is true in certain cases, in this article we have demon- 
strated that the finite temperature chiral phase transi- 
tion belongs to the 3d XY universality class even in the 
large N limit. In an earlier study we found similar results 
in two spatial dimensions [23j . In two dimensions contin- 
uous symmetries cannot break at any finite temperature 
However an 0(2) symmetry is special and a phase 
transition in the Berezinski-Kosterlitz-Thoulcss (BKT) 
universality class is possible [2^]. The BKT prediction 
for our model is that % ~ L 2- ^ where r\ is a function 
of temperature < r)(T) < 1/4 when T < T c . At the 
critical temperature r\ = 1/4. We found all this to be 
true even at large N. Interestingly, according to Witten 
a large N mean field theory would find r\ ~ 1/N [2fij |. 
While our results agree with this observation at a fixed 
temperature T, we find that if T/T c is held fixed then rj 
approaches a non-zero value in the predicted range show- 
ing that the mean field approach again breaks down in 
the large N theory close to the phase transition. 

Our study shows that the large N limit may not al- 
ways be able to suppress infrared fluctuations close to 
second order phase transitions. In a sense this is an in- 
dication that the perturbation expansion starting from a 
mean field solution breaks down. In other words a care- 
ful analysis of the FTMFT should be able to reveal this. 
This has not yet been done and is a useful project for 
the future. It can help us classify the types of large N 
models where mean field theory can break down. 
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